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Abstract. Results of first-principles calculations of the Fe/GaAs/Ag(001) epitaxial 
tunnel junctions reveal that hybridization of interface resonances formed at both 
interfaces can enhance the tunnelling anisotropic magnetoresistance (TAMR) of the 
■ systems. This mechanism is manifested by a non-monotonic dependence of the 

<— | , TAMR effect on the thickness of the tunnel barrier, with a maximum for intermediate 

thicknesses. A detailed scan of fey-resolved transmissions over the two-dimensional 
Brillouin zone proves an interplay between a few hybridization-induced hot spots and 
a contribution to the tunnelling from the vicinity of the T point. This interpretation 
is supported by calculated properties of a simple tight-binding model of the junction 
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, which reproduce qualitatively most of the features of the first-principles theory. 
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PACS numbers: 73.40.Sx, 75.70.Tj, 85.75.Mm 
1. Introduction 

Systems and devices with applicability in spintronics include traditional magnetic 
multilayers and tunnel junctions consisting of several magnetic layers [U [21 [3] as well 
as more recent artificial structures containing only a single magnetic part [HE]. The 
former exhibit the well-known giant and tunnelling magnetoresistance effects arising 
due to changes in the mutual orientation of magnetization directions, the latter are 
featured especially by properties driven by the spin-orbit (SO) interaction. The 
so-called tunnelling anisotropic magnetoresistance (TAMR) was observed in systems 
FM/I/NM, where the FM denotes a ferromagnetic metallic electrode, the NM denotes 
a non-magnetic metallic electrode, and the I denotes a non-magnetic insulating 
(semiconducting) barrier [6], [7] . 

The TAMR effect observed in conventional Fe/GaAs/Au tunnel junctions [7J 
proved to be rather weak which prevents its direct use. Similarly, a complementary 
phenomenon, namely, the current (or voltage) induced spin-transfer torques [4], which 
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might be employed for magnetization switching in the FM/I/NM devices, requires 
a large TAMR effect as a necessary prerequisite for efficient spintronic devices [8]. 
Recent attempts to enhance the TAMR values included, e.g., modification of the surface 
structure of the semiconductor layer during the growth process [9] or the use of an 
antiferromagnetic metal instead of the ferromagnetic one in the magnetic electrode [5]. 
The latter approach yields strongly enhanced TAMR values, which, however, could be 
observed only at low temperatures, despite the much higher Neel temperature of the 
antiferromagnetic material. 

The TAMR attracted interest also on theoretical side. A number of various 
topics were addressed in the framework of phenomenological models. These approaches 
discussed the role of the anisotropic density of states of the FM electrode [6], the 
symmetry properties and the interplay of the Rashba and Dresselhaus contributions to 
the SO interaction, the effect of external magnetic fields and applied bias voltages [7] [10] . 
etc. Ab initio calculations of the TAMR were carried out for tunnel junctions with Fe and 
Cu electrodes separated by vacuum [IT] and GaAs p2] as well as for the Fe/GaAs/Au 
trilayers [13]. The first-principles studies proved that interface states (resonances) 
formed at the FM/I interface play an important role in the TAMR phenomenon. 

All existing theoretical studies ascribe the TAMR primarily to electronic properties 
of the FM electrode, the tunnel barrier and their interface, whereas the NM electrode and 
the I/NM interface are considered to be of secondary importance. In the present paper, 
we show by means of first-principles calculations for the Fe/GaAs/Ag system and by 
the study of a simple tight-binding (TB) model Hamiltonian, that the standard picture 
of the TAMR is not generally valid. In particular, we find that interface resonances at 
the I/NM interface-when hybridized with those at the FM/I interface-can yield high 
TAMR values. Moreover, we predict a non-monotonic dependence of the TAMR value 
on the thickness of the tunnelling barrier and discuss its physical origin. 

The paper is organized as follows. The structural model of the Fe/GaAs/Ag 
junction and the ab initio techniques employed are summarized in section [2] while the 
obtained results and their discussion are presented in section El the properties of the 
GaAs/Ag interface are presented in section [37T1 and those of the Fe/GaAs/Ag system 
are contained in section 13.21 The TB model is formulated and studied in section H] and 
the summary of the main results is given in the last section. 

2. Model and methods 

The structural model of the Fe/GaAs/Ag(001) tunnel junction represents a simple 
generalization of the model for the Fe/GaAs/Cu trilayers [12]. The slab of the zinc- 
blende (zb) structure of GaAs is attached epitaxially to semiinfinite leads of the body- 
centered cubic (bcc) Fe and the face-centered cubic (fee) Ag; the atomic planes are 
the (001) planes of all three structures. The structure of the Fe/GaAs interface has 
been assumed without any reconstructions and layer relaxations on the basis of the 
ideal ratio of the lattice parameters a z \,/a\, cc = 2, that is satisfied with a good accuracy 
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for GaAs and Fe. The structure of the GaAs/Ag interface was also assumed without 
reconstructions; it employed the ideal ratio a z ^/a{ cc = a/2 and the mutual rotation of 
both bulk structures by tt/4 around the common (001) axis. The interplanar distance 
between the adjacent atomic planes of the GaAs and Ag parts was set equal to the 
arithmetic average of the distances in the bulk GaAs (a z b/4) and the bulk Ag (af cc /2). 
The z axis of the coordinate system is perpendicular to the atomic planes, while the x 
and y axes coincide respectively with the [100] and [010] directions of the fee lattice, 
i.e., they point along the [1, ±1, 0] directions of the bec and zb lattices. The GaAs slab 
contains n As-atomic (001) planes with As-termination on both sides; the investigated 
systems are thus abbreviated as Fe/As(GaAs)„_i/Ag(001). 

The electronic structure of the systems was calculated by self-consistent scalar 
relativistic tight-binding linear muffin-tin orbital (TB-LMTO) method [HI [15j [16] . The 
calculations were performed using the local spin-density approximation with the Vosko- 
Wilk-Nusair parametrization of the exchange-correlation potential [17]. The LMTO 
valence basis Ls, where L = (£, m) is the orbital index and s =f, 4- is t ne s pi n index, was 
limited to i < 2 and the so-called empty spheres were used for an efficient treatment 
of the open zb structure. The low-lying Ga-3d orbitals were treated as core orbitals 
so that the valence basis comprised the Ga-4c? orbitals [18] . This choice leads to a 
good description of the bulk bandstructure of the GaAs concerning both the band gap 
(around 1.2 eV) and the valence bandwidth (around 6.8 eV), which agree quite well 
with measured data [T9] . 

The study of the TAMR was based on the conductances evaluated in the current- 
perpendicular-to-the-planes (CPP) geometry. For this purpose, the scalar relativistic 
TB-LMTO Hamiltonian was completed by adding an on-site SO term having a simple 
£L- S form, where the atomic-like SO parameters £,r£,ss> of the atom at the lattice site R 
were calculated from the self-consistent electronic structure [20] . The matrix elements 
of the SO term have the form 

#fL,K'LV = S R R>hL<Bl,sA Ls \ L ■ U + SU\L' S '), (1) 

where the symbols L and S denote respectively the orbital and spin angular momentum 
operators. The direction of the FM electrode magnetization is given by two angles 9 
and 0, or by the unit vector n = (sin 9 cos </>, sin 9 sin 0, cos 9). These quantities enter 
the Hamiltonian via the operator U in H so (CQ), which acts only on the spin indices 
s,s' and which is represented by the spin- 1/2 rotation matrix U = D^ 1 ' 2 ' ((f) ,9,0) = 
exp(— i0cr 2 /2) exp(— i9a y /2), where the a y and o~ z are two of the Pauli spin matrices [21]. 
Note that the form of ([I]) corresponds to the global rotation of the spin quantization 
axis, whereas the coordinate system for the orbital motion remains unchanged. 

The CPP conductances were evaluated in the TB-LMTO Kubo-Landauer formalism 
[22l [23] as averages of the fey-dependent transmissions T(kn) over the two-dimensional 
(2D) Brillouin zone (BZ). The mixing of both spin channels due to the SO interaction 
has been implemented similarly to the case of non-collinear layered spin structures 
[22| [23]. For thicknesses of the GaAs barriers relevant in experiments (n > 20), the 
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Figure 1. The Bloch spectral functions (BSF) of the system Ag/(GaAs)is/Ag(001) 
for fen = f : (a) The layer-resolved BSF for the energy E located 1.5 mRy below the 
Fermi level. The dashed vertical lines denote the left (Ag/Ga) and right (As/Ag) 
interfaces, (b) The BSF's of the As interface atom as functions of energy, resolved 
according to the irreducible representations of the point group C% v . The vertical line 
denotes the position of the Fermi level of fee Ag. 



tunnelling current is carried mainly by the states with fey-vectors from a small central 
region of the whole 2D BZ. A sufficiently dense mesh of sampling points (corresponding 
to 5.4 x 10 5 fey-vectors in the full 2D BZ) has been used to get reliable values of the 
CPP conductances. 

3. Results and discussion 

3.1. Electronic structure of Ag/GaAs/ Ag(00l) systems 

In recent theoretical studies, the Fe/GaAs(001) system has been attached to a 
hypothetical bcc(001) Cu electrode since the bulk bec Cu has a free-electron-like 
bandstructure and the GaAs/Cu interface has a featureless transmission function 
[TT1 li~2] . This setup is advantageous for investigations of the role of the Fe/GaAs(001) 
interface state (resonance) lying in the minority-spin channel at the Fermi energy. This 
interface resonance gives rise, e.g., to a reversal of spin polarization of the tunnelling 
current with applied voltage p2] and it can lead to a pronounced TAMR effect [TT] . 
Motivated by these facts, we have focused on electronic properties of the non-magnetic 
GaAs/Ag(001) interface prior to the study of the Fe/GaAs/Ag junctions. 

We have thus studied systems Ag/GaAs/Ag(001) with different thicknesses and 
both terminations (Ga or As) of the GaAs barrier. These studies were done first without 
the SO interaction. Since the tunnelling current is mostly carried by states with fen- 
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Figure 2. The fc||-resolved transmission of the Ag/As(GaAs)i4/Ag(001) system at 
the Fermi energy: (a) without the SO interaction and (b) with the SO interaction. 
The coloured scale of the T(k\\) for both panels is shown on the right. 



vectors in the vicinity of the f point, we paid special attention to k\\ = f. Figure [D^a) 
displays the layer-resolved Bloch spectral function of the system Ag/(GaAs)i 5 /Ag(001) 
for an energy slightly below the Fermi energy (E — Ep = — 1.5 mRy). One can see a 
clear indication of an interface state at the GaAs/Ag interface, i.e., at the As-terminated 
GaAs barrier. The amplitude of the interface state is maximal in the As layer adjacent to 
the Ag electrode. A similar interface state was found in the Au/GaAs/Au(001) system, 
located however 10 mRy below the Fermi energy, whereas in the Cu/GaAs/Cu(001) 
with bcc Cu, no such state appears, in agreement with previous studies p2]. It should 
be noted that no similar interface state was found at the Ga-terminated boundary of 
the GaAs barrier in a wide energy interval around the Fermi energy (inside the band 
gap of GaAs), irrespective of the electrode metal (Cu, Ag, Au). 

The origin of the interface state can be understood from the Bloch spectral functions 
of the boundary As site resolved with respect to the symmetry given by the point group 
of the interface, namely, the C*2v group. This group has four one-dimensional irreducible 
representations: Ai, A 2 , Bi, and B 2 [25], of which only the Ai (subspace spanned 
by orbitals s, p z , d z 2 and d x 2_ y 2) is compatible with the symmetry of propagating 
states of the Ag(001) electrode at the Fermi level. As can be seen in figure H^b), 
the interface state is entirely of the symmetry Bi (subspace spanned by orbitals p y and 
d yz ), which is incompatible with the propagating states available in the Ag electrode. 
This incompatibility is an important factor, since the semiinfinite metallic electrode acts 
essentially like a vacuum half-space in the formation of the interface state. 

In the Landauer picture of the ballistic transport, interface states do not contribute 
directly to the system conductance, since the latter is given solely by the transmission 
coefficients between the propagating channels of the leads [26] . However, if an interface 
state is coupled, e.g., by a weak interaction to the propagating states, it can become 
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Figure 3. The density of states in two neighbouring Ga- and As- atomic planes in 

the middle of the Fe/As(GaAs)io/Ag(001) system for the majority spin ( , blue) 

and the minority spin ( , red). The vertical line denotes the position of the Fermi 

level. 



a resonance with a non-negligible effect on the conductance. In the present case, 
the SO interaction provides such a coupling of the Bi-like interface state to the Ai- 
like propagating state, which follows from an analysis of the double group C 2v and 
its irreducible representations: all spin-orbitals belong to the single additional two- 
dimensional representation of this double group [27]. The influence of the interface 
resonances on tunnelling is especially strong in symmetric junctions with identical 
electrodes, where both resonances become hybridized across the tunnel barrier which 
enhances the conductance appreciably [28]. This phenomenon is illustrated in figure EJ 
where the fey-resolved transmissions T(fey) are compared for the symmetric junction 
Ag/As(GaAs)i4/Ag(001) treated without and with SO coupling. The pronounced 
enhancement of the T(fcy) in vicinity of the F point is clearly visible; the total 
conductance of the junction increases by one order of magnitude due to the SO 
interaction. This result indicates importance of the GaAs/Ag interface for the transport 
behaviour of the Fe/GaAs/Ag system. 

3.2. Conductances and TAMR of Fe/GaAs/Ag junctions 

Figure [3] shows the local density of states (DOS) of the central GaAs layer (two 
neighbouring atomic planes) inside the Fe/As(GaAs)io/Ag(001) junction. The shape 
of the DOS is bulk-like, with the band gap clearly formed around the Fermi energy 
and with negligible spin polarization. These features prove that the junction is in a 
tunnelling regime with metal-induced gap states significantly suppressed for this and 
higher barrier thicknesses. The tunnelling regime is also manifested by an exponential 
decay of the conductance with the increasing GaAs thickness n, plotted in figure BJa) 



Tunnelling anisotropic magnetoresistance of Fe/GaAs/Ag 



-4 



o 

CD 
O 



-8 



(a) 




M 


|x X 


-- M 


|y \, 


M 


|z \ 



40 



cc 
< 



20 



(b) 




20 



30 



40 



20 



30 



40 



Figure 4. Dependence of the transport properties of the tunnel junctions 
Fe/As(GaAs)„_i/Ag(001) on the thickness n: (a) the conductances for the Fe 
magnetization pointing along the x, y and z axis, and (b) the corresponding in-plane 
TAMR. 



for three orientations of the iron magnetic moment, i.e., along the x, y, and z axis. 
For a given n, the conductance is obviously sensitive to the magnetization direction 
which leads to the TAMR effect. The orientational dependence was studied in detail for 
magnetization directions in the atomic planes, i.e., as functions of the angle <p with the 
fixed value of 9 = tt/2. The resulting angular dependences, shown in figure [5] for n = 23 
and n = 33, reflect the two-fold rotation symmetry (point group C 2v ) of the system, 
in full agreement with previous calculations [13J and measurements [7] performed for 
similar Fe/GaAs/Au(001) junctions. 

The in-plane TAMR, defined from the angular dependence of the conductance 
C(8, 0) as TAMR = [<7(tt/2, tt/2) - C(tt/2, 0)]/C(tt/2, 0), is presented in figure S^b) as a 
function of the barrier thickness n. The calculated TAMR effect is quite large, exceeding 
10%, in reasonable agreement with the values calculated for the Fe/GaAs/Au system 
[13], but about two orders of magnitude stronger than the experimentally observed 
TAMR values [7J. Moreover, the calculated thickness dependence is non-monotonic 
with a maximum obtained around n = 27 corresponding to the GaAs thickness of about 
7.5 nm. In order to identify possible reasons for the calculated high TAMR values and 
the non-monotonic thickness dependence, additional analysis is needed. 

In general, discrepancy between the calculated and measured transport properties 
of epitaxial magnetic multilayers can often be ascribed to imperfect atomic structure at 
the interfaces. The correct treatment of structure defects on an ab initio level employs 
either supercell techniques [291 E0] or effective medium approaches [23J combined with 
a particular microscopic model of the structure imperfection. In order to get a rough 
insight into the sensitivity of the TAMR to the quality of interfaces, we adopted here 
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Figure 5. The conductance C(tt/2, </>) of the systems Fe/As(GaAs)„_i/Ag(001) as a 
function of the angle <f>: (a) for n — 23 and (b) for n = 33. 

Table 1. Dependence of the in-plane TAMR on the position of two atomic planes 
with disorder in the Fe/As(GaAs)26/Ag(001) junction. The first row corresponds to 
the ideal system. 

position TAMR (%) 

41.2 

inside Fe 39.4 

at Fe/As 3.7 

inside GaAs 41.4 

at As/Ag 4.7 

inside Ag 32.2 

at Fe/As and As/Ag 3.6 



a simplified approach. We have simulated chemical disorder in the system by a finite 
imaginary part e > of the energy arguments z = Ep ± ie of the potential functions in 
the TB-LMTO Kubo-Landauer formalism [22l [23] . This modification was used only in 
a few selected atomic planes of the whole system: in two neighbouring planes located 
at a single interface (Fe/As or As/Ag), at both interfaces, inside the GaAs barrier, or 
inside the metallic electrodes. The value of e = 5 mRy was used in all cases. The 
results for n = 27 (the thickness corresponding to the maximum TAMR effect in the 
perfect junctions) are collected in table [D One can see that the disorder deep inside 
each part (Fe, GaAs, Ag) of the junction has only a minor influence on the resulting 
TAMR. However, the interface disorder reduces the TAMR value quite significantly, 
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Figure 6. The fc||-resolved transmissions of the Fe/As(GaAs)„-i/Ag(001) systems at 
the Fermi energy for 6 = tt/2: (a) n = 19, (f> = 0, (b) n = 19, = tt/2, (c) n = 27, 
= 0, (d) n = 27, = tt/2, (e) n = 47, = 0, and (f) n = 47, = tt/2. The coloured 
scales of the T(kii), shown on the right, refer to both plots with the same n. 



which proves that at least a part of the difference between the large theoretical values 
of the TAMR and the much weaker observed effect is due to the interface roughness. 
Moreover, both interfaces influence the TAMR effect to a similar extent, see table [TJ 
which indicates that they are of equal importance for the calculated trend of the TAMR 
(figure 1Kb)). 
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Further insight into the obtained results follows from the fey-resolved transmissions 
T(kn). Figure |6] shows this quantity for three barrier thicknesses, n = 19, n = 27, and 
n = 47, and for the magnetization directions along the x and y axis. Only a small 
region around the T point is included in the figure while the rest of the entire 2D BZ 
(defined by ab cc |^,j/| < tt/v^ ~ 2.22) is unimportant for the tunnelling. The maximum 
TAMR {n = 27) corresponds to a few hot spots in the T(k\\) plots, see figure Etc, d); 
their positions and contributions to the total conductance are strongly sensitive to the 
direction of the iron magnetic moment. The lower TAMR values for thinner barriers 
are due to non-negligible contributions of bigger 2D regions to the total conductances, 
see the case of n = 19 in figure El^a, b). On the other hand, the lower TAMR values for 
very thick barriers (n > 40) are due to the same hot spots with non-zero fey-vectors as 
found for intermediate thicknesses (n m 27) accompanied by another pronounced local 
maximum of the T(fey) in the T point, see the case of n = 47 in figure EJ^e, f). The 
contribution of the latter maximum to the total conductance is little sensitive to the 
magnetization orientation which explains the reduction of the TAMR effect. 

The presence of the hot spots in the fey-resolved transmissions is undoubtedly an 
important factor contributing to the high TAMR and to its non-monotonic dependence 
on the GaAs thickness. Similar hot spots appeared in various magnetic tunnel junctions 
[28| |3T| |32| [33] and were shown to be a direct consequence of a hybridization of two 
interface resonances across the barrier [28]. In the present case, an interplay of these 
hybridization-induced hot spots with the contribution of the F point represents a new 
situation, relevant especially for large thicknesses of barriers with a direct band gap, 
such as MgO or GaAs |3~T] . 

4. Hybridized interface resonances in a tight-binding model 

In order to assess the role of the presence and hybridization of the interface resonances 
on the TAMR, we have formulated a simple TB model and investigated its properties. 
The atoms are placed in positions of a simple cubic lattice with the lattice parameter 
a; the atomic planes are the (001) planes. The active part of the FM/I/NM junction 
comprises N + 2 atomic planes labelled by an index p, p = 0, 1, . . . , N + 1. The plane 
p = corresponds to a FM layer, the plane p = N + 1 corresponds to a NM layer, 
and the tunnel barrier is represented by the planes p = 1, 2, . . . , N. We assume a single 
orbital per site and spin and a spin-independent nearest-neighbour hopping between 
the orbitals. The hopping elements are different for pairs of atoms in neighbouring 
atomic planes (hopping t) and inside the same atomic plane (hopping t). The atomic 
energy levels of the tunnel barrier and of the NM layer are spin-independent while those 
of the FM layer are exchange- split. The direction of the exchange splitting is given 
by an in-plane unit vector n = (cos 0, sin 0, 0). The FM layer is also influenced by a 
Rashba-type SO splitting derived from the canonical form H so ~ (p x er) • u, where the 
P — (Px,P y ,Pz) denotes the momentum operator, the a = (a x , a y , a z ) are the Pauli spin 
matrices and the unit vector v = (0, 0, 1) is normal to the atomic planes. We assume 
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full 2D translational symmetry so that the Hamiltonian of the system can be written 
after the 2D lattice Fourier transformation as: 

H pSjP , a i(k\\) = 5 pp/ yif) - 2t[cos(ak x ) + cos(ak y ))6 ss >j 

— 6\p- p >ii$ss't + SppSp'pHgg, (fe||), (2) 

where the s and s' are spin indices, s,s' =t)i? which refer to the global (fixed) 
spin quantization axis along the (001) direction. The term nfj, comprises all on-site 
interactions in the p-th atomic plane and the last term describes the Rashba-type SO 
interaction in the FM layer (p = 0). For the NM layer (p = N + 1), the on-site term 
is given by h s Z +1 ^ = eNM^ss', and for the barrier layers (p = 1,2, . . . , N), it is given 
similarly as h^J, = tE^ss', where the parameters £nm and ee denote, respectively, the 
energy levels of the NM and barrier layers. The explicit form of the FM on-site term 
(p = 0) is given by 

'4°! = *S2 = ^ *S = "IT = ^ «*<-*). (3) 

where the parameters e-j. and denote the exchange-split energy levels. The last term 
in (E5 is given by 

i*$(fc B ) = flj5(fc,i) =o, 

fltJ(Jfe||) = #J°(fc||)* = - a[sin(aifc B ) + isin(aJfe !B )], (4) 

where the parameter a scales the Rashba-like SO interaction. The effect of the 
semiinfinite FM and NM leads has been simplified by adding energy- and fey-independent 
selfenergy operators to the on-site interactions of the NM and FM layers. These 
(retarded) selfenergies are given by 

V (NM) 1 r V^( FM ) V( FM ) 1 ( I \ 

Ks> = -^InmSss', S| t = E^ 4 = --( 7t + 7 J, 

S t™ } = _ ^(^t - li) ex P(-i0)> s | F t M) = ~^(7t - 7^) exp(i0), (5) 

where the parameters 7nm, 7t and 74, define the widths of the respective energy levels 
(resonance widths in the local spin reference system). 

The simplicity of the model allows one to choose easily its parameters in order 
to achieve the presence or absence of the resonance at either interface and for each 
spin channel. Here we confine ourselves to the case of one spin channel (spin- J,) out 
of the resonance and with no propagating states in the FM lead; the latter condition 
is obtained by setting 7^ = 0. We have considered four cases of the model. The first 
case, denoted 2, corresponds to the presence of two resonances: one at the FM/I 

interface in the spin-'f channel, the other at the I/NM interface. The particular values 
of the model parameters are: t = 0.48, t = 0.03, = 0.36, = —0.1, ee = 1.1, 
£nm = 0.4, 7| = 0.009, 7nm = 0.007, and a = 0.03. Note that a small asymmetry 
has been intentionally introduced in the parameters e^/eNM and 7t/7NM in order to 
simulate different properties of the FM and NM electrodes. The Fermi energy is set to 
zero, Ep = 0, which is located slightly below the bottom of the spectrum of the tunnel 
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Figure 7. Dependence of the conductance of the FM/I/NM model on the angle <j> for 
the barrier thickness N = 30 and for the case 2. 
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Figure 8. Dependence of the transport properties of the FM/I/NM model on the 
barrier thickness N for the four cases (for details, see text): (a) the conductance for 
= 0, (b) the in-plane TAMR. 

barrier, Vb = £b — 2t — 4t = 0.02. The second case, denoted as case IF, corresponds 
to the case of one resonance, located at the FM/I interface for the spin-t channel; its 
parameters coincide with the case 2 apart from the value of €nm = 0.6. The third 
case, denoted IN, describes the situation with one resonance, present at the 

I/NM interface. This case is obtained from the case 2 by setting the value of = 0.6. 
The last case, denoted as case 0, refers to the absence of any resonance; its parameters 
are obtained from the case 2 by changing its two parameters, namely, e-f = 0.6 and 
e N M = 0.64. 

The angular dependence of the conductance C(<f>) is plotted in figure [7]for the case 2 
and the barrier thickness N = 30. It is seen that the calculated dependence reflects the 
fourfold symmetry of the system (point group C*4 V ). The in-plane TAMR is thus defined 
as the ratio [C(7r/4) — C(0)]/C(0). The dependence of the conductances on the barrier 
thickness for all four cases and for <fi = is shown in figure E^a). One can see clearly the 
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effect of the resonances, pronounced especially for smaller thicknesses: the straight line 
for the case is slightly modified by the presence of a single resonance (cases IF and 
IN), whereas the hybridization of two resonances is manifested by a strong enhancement 
of the conductance (case 2). The corresponding thickness dependences of the TAMR 
are depicted in figure ffljo). It is seen that a sizeable TAMR effect is obtained only for 
the case 2 while the single resonances (on either side of the barrier) lead essentially to 
negligible TAMR values (cases IF and IN), similarly to the case 0. The hybridized 
interface resonances yield also a non-trivial dependence of the TAMR on the barrier 
thickness: small initial values for N < 10 are followed by a steep increase to a broad 
maximum for 20 < iV < 30 which is replaced by a final decrease for N > 40. This trend 
is qualitatively similar to that obtained for the Fe/GaAs/Ag system, see figure H^b). 

The different regimes of the thickness dependence of the TAMR can be related to the 
corresponding fey-resolved transmissions shown in figure [9]for the case 2 and three values 
of N. For N = 10, the total conductances arise from contributions of substantial parts of 
the whole 2D BZ which leads to a modest TAMR effect for thin tunnelling barriers. For 
iV = 20, the dominating contribution to the tunnelling is due to a narrow region (a hot 
spot), the position of which depends on the angle <fi. The sensitivity of this sharp single 
local maximum to the angle gives rise to enhanced TAMR values for intermediate 
barrier thicknesses. For iV = 40, the hot spots survive but are accompanied by a 
pronounced peak in the very centre of the 2D BZ, which is reflected by reduced TAMR 
values for very large insulator thicknesses N. This reduction is a simple consequence 
of vanishing of the Rashba term iff s ?(fey) in the limit of k x — > 0, k y — > 0, see (jlj). 
The obtained changes in the fey-resolved transmissions are in close analogy to the first- 
principles results, which corroborates the conclusions drawn in section 13.21 

5. Conclusions 

The calculated ab initio results for the Fe/GaAs/Ag(001) system with perfect epitaxial 
interfaces and the properties of a simple tight-binding model of the FM/I/NM junctions 
demonstrate that hybridized interface resonances can strongly influence the TAMR. In 
particular, they can lead to sizable TAMR values, especially for intermediate thicknesses 
of the tunnel barriers. The hybridized interface resonances can be thus added to 
the list of existing origins of the TAMR: the anisotropic density of states of the FM 
electrode [6], the interference effects of the Rashba and Dresselhaus contributions to 
the SO interaction [3, HD], and the interface states at the FM/I interface [11] . This new 
mechanism - if realized experimentally in a special junction - might also be employed to 
enhance the TAMR effect for applications in spintronics. Further open problems related 
to the presented results, such as, e.g., the effect of a finite bias, external magnetic fields 
or elevated temperatures, remain a task for future studies. 



Tunnelling anisotropic magnetoresistance of Fe/GaAs/Ag 



14 



III 

- (a) 

/ 










I,, 

- (b) 


I, 





) 





,iiii 
" (c) 

- 


1 

- 


1 


- 



iii 

" (e) 


III 


c 

1,1, 





,1111 

- (d) 


1 


i 


1 


,1111 
" (0 

c 

1 


I 

J 

I 



u 



0.86 



0.43 



i 



2x10"' 



1x10"' 



5.2x10 



2.6x10 



10 



10 



-0.5 





ak„ / ji 



0.5 



-0.5 





ak Y / k 



0.5 



Figure 9. The fc||-resolved transmissions of the FM/I/NM model in the case 2: (a) 
N = 10, cj) = 0, (b) N = 10, <f> = tt/4, (c) N = 20, = 0, (d) N = 20, = tt/4, (e) 
N = 40, 4> = 0, and (f) AT = 40, (j> = tt/4. The coloured scales of the T(k\\), shown on 
the right, refer to both plots with the same N. 
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